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We investigate the effects of inhomogeneous scalar field configurations on the elec- 
troweak phase transition. For this purpose we calculate the leading perturbative 
correction to the wave function correction term Z(ip,T), i.e., the kinetic term in 
the effective action, for the electroweak Standard Model at finite temperature and 
the top quark self-mass. Our finding for the fermionic contribution to Z(ip,T) is 
infra-red finite and disagrees with other recent results. In general, neither the order 
of the phase transition nor the temperature at which it occurs change, once Z(<p, T) 
is included. But a non-vanishing, positive (negative) Z(<p,T) enhances (decreases) 
the critical droplet surface tension and the strength of the phase transition. We 
find that in the range of parameter space, which allows for a first-order phase tran- 
sition, the wave function correction term is negative — indicating a weaker phase 
transition — and especially for small field values so large that perturbation theory 
becomes unreliable. 
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I Introduction 



The understanding of the electroweak phase transition has matured rapidly during recent years. 
The original work on restoration of gauge symmetries at high temperatures [|IJ and systematic, 
partial summation of perturbation theory |2j has attracted a lot of attention in the framework of 
the Standard Model, since the suggestion of a mechanism for baryo-genesis at the electroweak 
scale Already in the early eighties it became clear that non-abelian gauge theories contain 
massless degrees of freedom in perturbation theory at high temperatures Therefore, infrared 
problems complicate the straightforward setup of perturbation theory and a considerable effort 
was made recently to resolve them P-|jT7|. Up to Higgs masses of about 80 GeV, these difficul- 
ties may be cured by introducing improved propagators (c.f. e.g., |14|-[I7|). In the improved 
perturbation theory self-energy and self-mass corrections are included in the action from the 
beginning. They are determined by solving appropriate special cases of the Dyson-Schwinger 
equations, the gap equations 0, [16], [L7]]. Keeping so-called "hard thermal loops" and summing 



over "soft thermal loops" the leading infrared singularities can be shown to cancel out |T8 



The order of the phase transition depends crucially on the value of the magnetic mass of 
the gauge bosons, whose calculation requires non-perturbative techniques. Although different 
approaches suggest the same order of magnitude fl6], the determination of its value to 
within more then 10% accuracy would be desirable. As long as the magnetic mass is smaller 
than 0.07 T, where T denotes the temperature of the heat bath, the phase transition appears 
to be of first order in the range of Higgs masses accessible to perturbation theory. However, it 
is still an open question, whether the phase transition within the electroweak standard model 
is strong enough to account for the observed baryon asymmetry of the universe |20| . Due to its 
weakness the phase transition may even proceed via formation of sub-critical droplets |2lJ . 

All the above results have been extracted purely from the effective potential V e s((p,T), 
i.e., the effective action for vanishing derivative terms d^ip = of the scalar field if, which 
plays the role of the order parameter. In this paper we investigate whether the full effective 
action is indeed dominated by homogeneous field configurations and whether it is justified to 
neglect quantum corrections which give rise to derivative terms in the effective action. We 
wish, however, to point out from the beginning that, as a matter of fact, a consistent one-loop 
treatment of the phase transition requires the inclusion of the wave function correction term 
whose leading contribution occurs at this level in perturbation theory. Indeed, it was argued in 
refs. [22|— [25 that depending on the theory under consideration the expansion of the effective 
action around d^ip = might break down and non-perturbative effects become important. In 
addition we will see that quantities essential for the mechanism of critical droplet nucleation 
such as the droplet surface tension and the strength of the phase transition are modified by the 
presence of higher derivative terms. Although, in this paper, we will perform a perturbative 
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calculation there are techniques of averaging the action over a range of momenta without 



expanding it in terms of derivatives p6 



The outline of the paper is as follows: In section [TJ we will introduce some essential quantities 
for the description of first-order phase transitions, review their derivation for homogeneous 
scalar field configurations and extend the analysis to inhomogeneous field configurations. In 
particular, we explain the impact of the wave function correction factor Z(ip,T) on droplet 
nucleation rate and surface tension for a first-order phase transition. In section |T| we discuss 
the self-energy corrections for the various particles in the theory in the presence of a plasma 
which help to cut off infrared divergences. A general method to derive the kinetic term in 
the effective action developed in [27 is reviewed in section [IV] where we explicitly calculate 
the different contributions to the wave function correction term for the Standard Model. Our 
results are described in section [V], and our conclusions are presented in section [VI[ A discussion 
of the top quark self-mass and several useful integrals with their high-temperature expansions 
have been relegated to two appendices. 



II The wave function correction term 
A Decay of metastable states 

The entire dynamics of the phase transition is contained in the quantity 

Z fi [<p] = l[V<p][DW^}e W {S p [cp + <p,W^}} . (II.l) 

The path integral is performed over fluctuations <p around a classical field configuration <p(x), 
over vector and ghost fields, whose measure is collectively denoted by [IJWJ and over fermion 
fields with measure [Dip]. The exponent contains the classical action Sp at finite tempera- 
ture T = An additional term that vanishes for stationary field configurations has been 
neglected. 

Following the approach of [16], in a first step we integrate out all vector, ghost and fermion 
fields (but not the scalar field fluctuation (p) to arrive at an effective (coarse-grained) finite- 
temperature action Tp[(p] defined via 

Zf>\ip\ =: J [V0] exp {T p [ip + 0}} , (II.2) 
which may be expanded in powers of derivatives of tp: 

I>]= / \l(d<p?-V(<p,T)}+ [ fiz^T)^)". (II.3) 
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Each Z n may in addition be expanded in a power series in the coupling constants of the theory. 

The remaining integration over ip in ( fll.2|) will be carried out in the saddle point approx- 
imation. We collectively denote all possible Lorentz invariant derivative terms of (p with n 
derivatives^ by (dip) n and use the shorthand notation 

/. = / 

J/3 JO 



dr / d 3 x . (II.4) 



Throughout this paper we will use the imaginary time formalism [BR |2U[. Therefore boson and 
fermion fields satisfy periodic and anti-periodic boundary conditions in the imaginary time 
t = ixo, respectively. In momentum space this leads to integrals 

where u = 2ninT (u = (2n + VjitiT) for bosonic (fermionic) fields and n runs over all integers. 
The first term in (|II.3|) contains the effective potential V(<p,T) and the classical kinetic term, 
whereas the second summand incorporates derivative terms due to quantum corrections. 

For stationary fields, d T <p = 0, the effective action Tp\<p\ plays the role of the free energy 

F[<p,T\ = -±r fi [<p], (11.6) 

and (p is the order parameter of the phase transition. Usually first-order phase transitions are 
studied under the assumption that the effective action T^^] is dominated by homogeneous 
field configurations, i.e., the derivative terms proportional to Z n in (|II.3|) are neglected. Then 
the remaining path integral in ( 11.2 ) is carried out in a saddle point approximation around 
configurations Tp(x) which extremize the classical free energy 

F[tp, T]=J d 3 x \V{^ T) + l - |Vyf} . (II.7) 

Isotropic configurations Tp(r = \x\) obey the differential equation 

d^ + 2dJp_dV_ /-n_q / IIg x 
dr 2 r dr dtp 

which is obtained by varying ( [11. 7| ) w.r.t. (p. 

The generic potential for a first-order phase transition and for Tj, < T < T c has two local 
minima (c.f. Fig. [I]) one of which is metastable. Here T c denotes the critical temperature, 
where the minima are degenerate, and is the so-called barrier temperature at which the 

3 For simplicity we have not taken into account the fact that for given n different combinations of derivatives 
may have different prefactors Z n . 
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potential barrier between the minima vanishes. We choose the metastable minimum to occur 
at ip = and denote the position of the global one by <p — <p m m(T). The appropriate boundary 
conditions for the tunneling solution, which interpolates between the symmetric (<p = 0) and 
the broken (tp ^ 0) phase, are Tp'{r = 0) = and Tp{r — > oo) = [Q. In the thin wall 
approximation |3]J one finds 



1 — tanh 



V - r(tY 
d 



(II.9) 



where 

m = (h.10) 

is the droplet radius at temperature T. It depends on the surface tension a{T) which may be 
evaluated at the critical temperature T c 

Vmin(T c ) 



a = dcp^2V(cp,T c ) (11.11) 

and the potential difference between the two minima 

AV (T) = V(0, T) - V(^ min , T) . (11.12) 

The thickness d of the droplet wall depends on the detailed shape of the potential. For a 
polynomial of the form 

V(<p, T) = l -m\T)^ - ET(p 3 + (11.13) 

one obtains 

d=J\—. (11.14) 

By inserting ( [II. 9| ) into (|II.7| ) we find the free energy in the thin wall approximation as a sum 



of a surface and a volume term: 

4-ir Air 
F TW [^,T] = AnR 2 (T)a - y i? 3 (T) AV(T) = -oR 2 {T) . (11.15) 

Once the thin wall approximation breaks down, AV(T) becomes larger than the barrier 



between the minima, and QII.8D has to be solved numerically. Unfortunately, in the Standard 



Model the thin wall approximation is only marginally applicable for m# ~ 80 GeV 0. [16 



To second order in the saddle point approximation of (|II.2| ) we need to perform a Gaussian 
path integral over the fluctuations (p around the tunneling solution Tp. This correction appears 
as a pre-exponential factor A in the nucleation rate of critical droplets 



r(T) = AT 4 exp {-{3 (F[<p } T] - F[0, T})} 



(11.16) 
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It has been shown to be of order one for the electroweak phase transition || |TjJ. In the following 
we will therefore only be concerned with the saddle point configuration Tp(r) . 

The temperature T e which marks the end of the phase transition may be defined as the 
temperature for which the droplet nucleation rate becomes larger than the expansion rate of 
the universe, i.e., T(T e ) ^ H A {T e ) 1 where H{T) denotes the Hubble function. Using ([II.16|) this 
leads to the rough estimate (3 e F\(p, T e ] ^ 145. 



B Inhomogeneous field configurations 

Higher order derivative corrections will in general alter the results of the last section. Some 
recent publications [[22|, ^3] indicate that depending on the parameters of the theory under 
consideration they may even dominate the leading terms and the expansion breaks down alto- 
gether. But even in the domain where the perturbative results seem to signal that this is not 
the case, non-perturbative considerations yield results which are different from perturbative 
calculations for the effective potential alone |24|, |25||. In this paper we will restrict ourselves 
to the first term of the derivative expansion, i.e., Z(<p,T) := Zz((p,T), that corrects the scalar 
wave function. Here we will merely discuss possible physical consequences of a non-vanishing 
Z(ip,T) and leave the explicit calculations within the Standard Model to the forthcoming sec- 
tions. We would like to point out that a consistent one-loop analysis of the electroweak phase 
transition indeed requires the inclusion of Z(ip,T), since its leading perturbative contribution 
occurs at the one-loop level. 

To check whether homogeneous field configurations are important for the dynamics of the 
phase transition we wish to investigate how Z(ip,T) ^ affects the characteristic quantities of 
the phase transition. For this purpose we perform the scalar field transformation 

(p{r) = J dtpy/l + Z^T) (11.17) 

under which the free energy 

F[<p,T\ = Jd*xl [ V(<p,T) + 1 -(V& 2 } (11.18) 

has the same form as for Z(ip,T) = 0. The new potential V(ip,T) = V(ip(<p),T) is locally 
rescaled. Since = y/1 + Z > there is a one-to-one correspondence between (p and (p, i.e., 
( |1I.17| ) amounts to a local rescaling of the yj-axis in Fig. |l[ Minima and maxima still have the 
same potential energy V((p min ,T) = V(<p m i n ,T) and as a consequence the critical temperature 
T c , at which the two minima are degenerate, does not change either. Also neither the height 
of the barrier nor the amount of supercooling AV(T) = AV(T) change. Thus, once the thin 
wall approximation is valid for homogeneous field configurations it survives the incorporation 
of the wave function correction term. 
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However, the new surface tension 

a = j o d(p^2V{{p,T) = J o d<py/2[l + Z(<p,T)]V(<p,T) (11.19) 

may be substantially different from the surface tension a without wave function correction term. 
Hence, a possible measure for the effect of the wave function correction term on the dynamics 
of the phase transition is given by 

OzF := = — (11.20) 

F[<p, T] a 3 

where <p denotes the stationary isotropic solution which extremizes the corrected free energy 
( [II. 18 ). For a constant Z = 0.25, for instance, one would obtain the significant relative deviation 
of 5 Z F ~ 0.4. 

Another important consequence of a non-negligible Z-factor is the end temperature of the 
phase transition. Suppose that Z((p,T) > 0; this would result in an increased surface tension 
and the universe would have to supercool further to complete the phase transition: T e < T e . 
Therefore, the strength of the phase transition would increase, since ^ mi ^ Te ^ > gsaslZk) anc i 
could be more favorable for baryogenesis than anticipated without wave function corrections. A 
rough estimate based on the potential ( |11.13| ) and the assumption of an approximately constant 
Z yields 

^ min J Te) = |i + - [(i + zy - i] ) ^ min(Te) , (n.21) 

where we assumed that X ™Jr^ « 1- Eqn. ( |H.21|) shows that a positive wave function 
correction factor Z would enhance the strength of a first-order phase transition whereas a 
negative Z would decrease it. 



Ill Self— energy corrections and improved propagators 

Finite-temperature field theory is known to develop infrared singularities once massless de- 
grees of freedom are present. Since in the symmetric phase all masses are essentially zero, 
straightforward perturbation theory breaks down. This problem can be avoided by using an 
improved perturbation theory with propagators that include finite-temperature self-energy and 
self-mass corrections generically denoted by E n (x) in this paper. In general, the 



full propagator V n {k) is determined by the Dyson-Schwinger equation 

v-\k) = v~} ) (k) + ti: n (k), (in.i; 
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where n labels the different fields in the theory and T> n fi(k) is the tree-level propagator. The 
Dyson-Schwinger equations may be solved perturbatively. To the order to which we are calcu- 
lating, the vertex functions can be replaced by the tree level couplings gi and we are left with 
a one-loop diagram for E n . 

In the following we will simply state the known self-energies £„ = n s (fc) for scalars and 
E n = —U v (k) for vectors in the electroweak Standard Model. The top quark self-mass S n = 
S f (k) is evaluated in appendix [A] in the limit of small k. 

The starting point for our explicit analysis is the simplified SUl(2) Standard Model (SM) 
with vanishing hypercharge gauge coupling g'. In this approximation the SUi(2) gauge bosons 
are degenerate in mass and the Weinberg angle Qw is neglected. The corresponding Lagrangian 
is 

£ ^gauge ^Higgs ^fermion ^GF ^ghost (III. 2) 

with 

£ gaU ge + £ffig gs = ~\f; v F^ v + (D,$) t (D"$) - fj 2 - A ($t$) 2 , (III.3) 



— dp — ig-W« , (III.4) 



/r < and 

r" 
~2 

where the r a , a = 1,2,3, denote the three Pauli matrices. The field $ is an SUl{2) doublet 
parameterized by four real scalar fields: 



$(x) = —( Xi(^)+%(a;) | / m5 \ 
U y/2 { V (x) + h(x)+ixs(x) 1 ' 

where h is the Higgs field, Xa, a = 1,2, 3, are the three Goldstone bosons and ip is a real 
background field. 

In the fermionic part of the Lagrangian we neglect all lepton and quark Yukawa couplings 
compared to the top quark Yukawa coupling f t . Hence 

Acrmion = V^D^L + ft (t L , h) («T 2 $*) t R + h.C. (III.6) 

The gauge fixing and ghost Lagrangians are given by 

^GF — ~7^ J ~ a,J~ 

2 £ (III.7) 

■^ghost = C a Ai ab Cij 

with T a (W) and M. defined via 

8 



M ab = __jrb^ W u 

OUJn 



uj=0 



T a = d,WZ- l -g&Xa 



(III.8) 
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where W w denotes the result of an infinitesimal gauge transformation U (u>) = 1 — iT a u a on the 
gauge field W. Throughout this paper we will work in Landau gauge £ = 0. 

To one-loop the full scalar propagator can easily be read off from ( [111.1| ): 

*W*) = ~ 2 m n m (IIL9) 

with m^ 2 = X(3(f 2 — v 2 ), m^ 2 = X((f 2 — v 2 ) being the tree-level scalar masses. Due to 
the breakdown of Lorentz invariance in the presence of a heat bath the full vector propagator 
(mg? = g<p/2) 

= fo? P L 



m$ 2 -n T (k) T k 2 -Zm$ 2 -Tl G (k) G 



involves projections onto the direction of propagation (Pq V ), onto the component of the 
heat bath flow perpendicular to (-P^) and onto the remaining two directions (P^)- For 
explicit definitions and properties c.f. ||16|| . In terms of these projectors the self-energy tensor 

W v (k) = Tl L {k)P£ u + U T {k)Pf + Tl G {k)P£ u + n s {k)S^ (III.ll) 

also involves a traceless projector S which does not contribute to T>^ v in Landau gauge. Note 
that we keep the third term in (LLL.lUp although at first sight it seems to vanish for our gauge 



choice £ — > 0. Eventually one encounters singularities in the gauge fixing contribution of Z(ip, T) 
for £ — > which are precisely canceled by this term. 

In the limit k — > the self-energies give rise to the plasma mass corrections 

9 2 + \\ + \tf)T 2 

2 T 2 (IIL12) 

1 2 9 4 T 2 . 

The magnetic gauge boson plasma mass 5mw,T of order g 2 is a non-perturbative feature of the 
theory which was first predicted in [|J and derived in e.g. |Ll| [Tj]]. The factor 7 is expected to 
be of order one as confirmed by lattice simulations and other non-perturbative methods ]19| . 



1W0) 


= 6m l,x = 


(l6' ( 


Hx(0) 


= 5m 2 W)L = 


11 


n T (o) 


= bvr? WT = 


1 

9tt 2 



Though fermionic contributions to the wave function correction term are not expected to 
introduce new infrared singularities, the top quark self-mass S t is evaluated in appendix [A[ 
We do not find any plasma mass correction to first loop-order. This is however not surprising, 
since at least in the symmetric phase m t = is protected to any finite order in perturbation 
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theory by chiral invariance. Instead we find a chemical potential which is non-vanishing in 
both phases of the theory. We will however neglect this one-loop correction to the tree-level 
propagator. In the next section we will find, that the (unimproved) top quark contribution to 
Z(ip,T) is indeed infrared-finite and furthermore negligible compared to other contributions. 
Therefore there is no need to improve perturbation theory by including this plasma induced 
effective chemical potential into the top quark propagator. 

The improved propagators for the bosonic degrees of freedom are now given in terms of the 
full masses 

m \ = m f )2 + 5m 2 . (III. 13) 

The tree-level connection between the couplings and zero-temperature masses is given by g = 
2m w (T = 0)/v, f t = V2m t (T = 0)/v and A = m 2 H (T = 0)/2v 2 where v = sJ-^/\ ^ 246 GeV. 



IV Local momentum expansion 



Previous work on terms quadratic in spatial derivatives in the effective action include |27| and 
p3|-||43||. We will follow the approach of Moss et. al. which is based on the local momentum 
space method of 



We start with a short review of the methods developed in [2?| and generalize them by 
including plasma masses. The calculation of the one-loop contribution to the wave function 
correction term in the effective action only requires the part of ( |111.2| ) which is bilinear in the 
quantum fluctuations^] h, Xi, W, c and t: 



Co 



U ab W;A^W b p + 5 ab c a A c c b + \h^h + W xAxX 3 + tA t t 



(iv.i; 



with 



AJx) 



A t (x) 



d p d p + m 2 w )g 



[IV 



1 - X - j d»d v 



-{d'd p + ^m 2 w + 2N^ 
id + m t 



(IV.2) 



and N„ = dn In ip. The one-loop corrections to the effective action may then be written as 



rW = -XXtr(lnA„(x)) 



(IV.3) 



Of course, h and \i are exactly the fluctuations which are collectively denoted by (p in (II. 1) 
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where n runs over all fields in the theory. As seen above, A n is an operator of the general 
form A n (x) = D n + M n (x) with D n being a differential operator and M n (x) a mass termQ 
The coefficients a n are fixed by the statistics and the number of degrees of freedom of the 
corresponding field. In particular, we have a v = a Xj = |, aw = § and a t — — 1. 

The effects of the plasma are taken into account by using the inverse propagators^] 



A n = (D n + M n + E n ) - E n = A; - E n {x) 



(IV.4) 



where S ra denotes all plasma corrections for the inverse propagator, and may in general depend 
on x via (p. Therefore, to one-loop, ( [IV .3|) is equivalent to the modified relation 



r (i) 



-a n tr 



(IV.5) 



The situation for gauge fields is slightly more complicated due to the existence of different 
plasma masses for transverse and longitudinal degrees of freedom. According to ( |111.10| ) , ( |i 1 1 . 1 1| ) 
the operator 



A w = (A' WfL - 5m 2 WL P L ) + (A' WiT - 5m 2 WT P T ) + A W , G + . . . , 
where the dots indicate terms depending on N^, splits up into three contributions: 



(IV.6) 



a; 



W,L 



A 



W,T 



(d p d p + 
(d p d p + 



m WL \ P L 



in 



Pt 



(IV.7) 



A 



G 



In Landau gauge the final result does not depend on the last term and thus we need no counter 
terms for Aq. 

The first step in evaluating ( |1V.5| ) is the definition of a (finite-temperature) Greens function 
G n (x, x') for each A' n via 

A' n G n (x,x') = 8 X0 , x ,5(x-x") . (IV.8) 

If we furthermore assume that M n may be written as M n (x) = m n + M n (x), where m n is 
independent of x and M n is independent of m n , we arrive at 



dm r 



-a n / tr 

Jp 



1 - S r 



d 



(IV.9) 



5 Note that M n (x) actually corresponds to a mass term only for fermions but to a mass squared term for 
gauge bosons and scalars. 



'One may equivalently use improved (inverse) propagators from the beginning in ( [V. l| ) . In this case one 



has to add appropriate plasma mass counter terms to (1V.1) in order to compensate for these corrections M 
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where the functional trace is accounted for by the integral. Hence the knowledge of the G n (x, x') 
will enable us to derive the one-loop contribution Y^> to the effective action via integration 
w.r.t. m n . 

Since we only need to know theTj G(x,x') in the limit x — > x', we expand the functions 
M'(x) := M(x) + into a Taylor series around x = x'\ 

OO 1 



with y = x — x' and 



M' (x 



1=1 «' 



<9'M'(x) 



Inserting this into ( |1V.8| ) and Fourier transforming w.r.t. x (denoted by a tilde) 

d 3 k 



(IV.11) 



aO = E -a I ^— 3 e^~^G(k, x') , (IV. 12) 

k P J ( 27r ) 



we arrive at 



[D + M'(x') + g T[ K,..,M') dKi _ dk J G(k,x') = 1 . (IV.13) 

Note that the sum in ( |IV. 12|) may be over bosonic (ko = 2ninT) or fermionic (ko = (2n + l)iTiT) 
Matsubara frequencies, depending on the type of inverse propagator A' used to define G(x, x'). 
We may now expand G(k, x') and M'(x') in a series w.r.t. the number of derivatives^: 

oo 

M'(x') = £i'M'« (ar 7 ) 



(IV.14) 



One immediately obtains 



and 



G(k,x f ) = Y,G (j \k,x') . 

3=0 



&°\k, x') = ^ (IV.15) 



G^(k,x') = -&°Hky)j2±jM'^(x') 8 \ GV~'\k,x') . (IV.16) 

This iterative solution of ( [IV. 8| ) may now be used to determine the Greens functions and 
subsequently the effective action via ([IV .9|) order by order in the derivative expansion ( |PV.14|) . 

7 For simplicity we will suppress the index n from now on. It is clear that the following steps have to be 
carried out for each degree of freedom labeled by n, separately. 

8 For scalar and fermion fields the derivative expansion of M' is trivial. Eqn. ( IV. 2 ) however shows that it is 
non-trivial for vector and ghost fields. 
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The &°\k,x') yield the effective potential whereas the & 2 \k,x') will give the desired kinetic 
term in the one-loop effective action: 

1 r d 3 k 



rW - _ V a f dm f lim V - / ^ J(k y -ky) tr 
i kin - l^anj rfm / lim> / e tr 



(IV.17) 



Although ( |1V.16| ) and ( |i V . 1 7| ) provide a well-defined prescription for calculating the wave 
function correction, it still requires rather lengthy calculations. We have used the symbolic 
manipulation package Mathematica to accelerate the computation. All that is needed are the 
(improved) zeroth order propagators G^'(k, x) as given in ( |111.9| ), ( [111. 10| ) as well as the standard 



tree level top quark propagator, and the space-time dependent pieces M'Jp{x) which may be 
read off directly from the results of section |TTT| . As explained in section |IJ we will only take 



into account stationary fields, i.e., d^ip = 0, and use standard integration techniques for the 
finite-temperature momentum integrals. All contributions split up into a T-independent and 
a T-dependent piece as is generally the case in thermal field theory. 

The temperature dependent scalar contribution turns out to be 

ZStta T) = + -> + ^ + . (IV.18) 

The corresponding T-independent part is UV-finite and given by 



Z^M- AV ( 3 + 35m ' + 1 + ^1 (TV 10) 

-j^[^ + -^r + ^2+^r) ■ (IV ' 19) 



Since it is subleading in the high-T expansion compared to Zj^^cp, T) we will neglect it in 
the following. The T-dependent gauge boson one-loop contribution turns out to be 



3o 2 T / m m 5m (0)2 5m (0)4 41m (0)2 1 
Z^Uv, T) = ^ - + - + — I ■ (IV.20) 

s "" r 4vr \32m| 96m| 16m 5 T 48m| m T 1 



The zero-temperature contribution is UV-divergent and has to be renormalized. It can be seen 
to be of the form |?7 

Z^M ~ ^ (^) , (IV.21) 



where m ren denotes the renormalization scale. Compared to (|1 V .2U|) this part of the wave 



function correction term is again subleading in the high-T expansion and may therefore be 
neglected. 

We now turn to the one-loop top quark contribution Z fermion (y?, T). In the high-temperature 
expansion its leading order T-dependent piece is given by: 

d°ion(^ T) = ^ (l + 3 lE + 3 In ^) , (IV.22) 
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which is logarithmically IR-divergent. However, in this case the renormalized-temperature 
contribution which could be neglected for the other fields comes to the rescue. It is given by 

ZlTXM = JL fl - 41n , (IV.23) 



327r 2 \ rn T 

and therefore cancels exactly the IR-divergence of the T-dependent part. The full fermionic 
contribution is of the form 

Zfermion^, T) = ^ (V + 12 7E + 12 In ^\ . (IV.24) 

We would like to point out that our results for the bosonic as well as the T-independent 



fermionic contributions to Z((p,T) agree with the findings of [27| for vanishing plasma masses 



However, the fermionic T-dependent piece ( IV.22j ) disagrees with that of |27j where a l/m ( 



infrared singularity (which can not be canceled by the T-independent part) was found. The 
reason for this discrepancy appears to be that in |27| the kinetic energy was calculated following 



( |IV.17[) by summing over bosonic frequencies ko instead of fermionic frequencies as required 
by the Fourier transformation ( |1V.12| ) of a fermionic Greens function. Our results show that 



all infrared singularities in the symmetric phase are removed by the plasma masses to leading 
order in the high-T and the loop-expansion. 

The full wave function correction term is now given by 

Z(ip, T) = Z vectoI (ip, T) + Z {cimion ((p, T) . (IV.25) 

Note that here Z sca i ar is not included. We wish to emphasize that it would be inconsistent to 
add the scalar contribution for the calculation of physical quantities like the free energy, the 
droplet tension or droplet nucleation rate, following the methods described in section 0.0. The 
scalar contributions to the effective action are fully taken into account by solving the bosonic 
path integral ( |11.2| ) in the saddle point approximation. However, for checking the range of 



validity of the derivative expansion ( fT^f ) one might as well integrate out the scalars in the 
same way as the other fields and consider the quantity 

Z(tp, T) = Z sca i ar (yj, T) + Z vectoi ((p, T) + Z ievmion (ip, T) . (IV.26) 



V Results 

A short look at the various contributions to Z(ip,T) shows that it is completely dominated 
by its gauge boson contribution Z vector ((p, T) for temperatures close to the critical temperature 
T c . This is demonstrated in Fig. 0, where Z and all its contributions are plotted for physically 
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realistic values of the parameters m#(T = 0), m t (T = 0) and 7 at T = T c . We also note that 
an inclusion of Scalar would be numerically insignificant. 

Since the perturbative expansion of the wave function correction term starts at the one- 
loop level, the size of Z((p,T) is a direct measure for the convergence of perturbation theory. 



Furthermore, for values of <p and T for which \Z(<p,T)\ ^ 1 or Z(ip,T) 



<: 1, we have to 



expect that other higher derivative terms of the effective action in (|11.3| ) could also become 
important and may no longer be neglected. It is already known [IB], that there is only a 
small window of parameters close to the current experimental bound on the Higgs mass, where 
perturbation theory can be trusted as far as a calculation of the effective potential is concerned. 
In particular, for Higgs masses above approximately 85 GeV higher order loop-contributions to 
the effective potential become dominant and perturbative statements on the phase transition 
are questionable. A first order phase transition can reliably be found for m H ^ 85 GeV and 
small 7. It is altogether excluded for 7 2 taking into account the current lower bound on the 
Higgs mass of m# ^ 63.5 GeV [[44 1. 

In Fig. we have plotted the curve which divides 7-</2-space into two regions. Above it 
\Z((f,T)\ < 1 and the wave function correction term is smaller than the classical kinetic term 
of the effective action. Below this curve \Z(ip,T)\ > 1, perturbation theory breaks down and 
higher derivative terms might dominate the dynamics of the phase transition. We have used 
\Z(<p,T)\ ~ |Z vect0 r(y2, T)\, which is a good approximation for all temperatures relevant for the 
phase transition. We find that for the whole range of magnetic masses which permit a first-order 
phase transition (taking into account the experimental lower bound on the Higgs mass), there is 
an interval of small ^-values for which the wave function correction dominates over the classical 
kinetic term and is furthermore negative. Therefore, in this range the field transformation 
( |11.17| ) yields an imaginary scalar field (p and a complex effective potential V(<p>,T). This in 
turn renders the methods of section || inapplicable. However, a rough estimate, using an 
averaged (constant) wave function correction factor < Z(ip,T) >= J ¥ ' min dip Z(tp,T), gives 
e.g. for 7 = 1, rriH = 85 GeV and rrit = 160 GeV the value < Z(<p,T) >~ —0.4. This in 
turn yields an uncertainty of the derivative expansion of approximately 50% (c.f. (|II.20|) ) and 
a weakening of the phase transition strength by 3%. 

At the one-loop level, we are thus left with the following situation: For small magnetic 
masses (7^2), i.e., for the range of 7 which, in combination with the current experimental 
lower bound on the Higgs-mass, allows for a first-order phase transition, neither perturbation 
theory nor the derivative expansion seem to be reliable for small scalar field values. In particular, 
for those values of (p which are important for the determination of the dynamics of the phase 
transition, < (p < (p m i n , the one-loop wave function correction term dominates over the 
classical kinetic term in the effective action. Therefore, in this range of parameter space, it is 
not possible to make physically reliable predictions at the one-loop level. In order to clarify the 
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dynamics of a first-order phase transition, clearly, the calculation of higher loop-corrections to 
Z((p,T) as well as higher derivative terms of the effective action are necessary However, since 
the one-loop result for Z((p, T) is negative for the whole range of < ip/T < 1, we expect that 
its inclusion will weaken the first-order electroweak phase transition, as explained in section [H[ 

For larger values of the magnetic mass, 7^2, the perturbative calculation of the wave 
function correction term appears to be reliable. On the other hand such magnetic masses 
correspond to a part of parameter space, where the phase transition is expected to be of second 
order [p^| , and baryo-genesis can not appear at the electroweak scale within the mechanism of 

B, & 

VI Conclusions 

We have investigated the impact of the scalar field wave function correction term Z((p,T) on 
the electroweak phase transition. Since Z(ip, T) is the quantum correction to the classical scalar 
field kinetic term in the effective action, its knowledge is crucial for an understanding of the 
influence of inhomogeneous scalar field configurations on the dynamics of this transition. 

It turned out that generally a positive Z-factor would increase the strength of a possible 
first-order transition whereas a negative Z would decrease it. A calculation of Z(<p,T) within 
the Standard Model to one-loop revealed that it is fully dominated by its gauge boson contribu- 
tion for temperatures close to the critical temperature T c and well above the barrier temperature 
TV Furthermore we find that the one-loop result for Z(<p,T) is negative for the whole range 
< < 1 and T ~ T c . However, for small magnetic gauge boson masses, 7^2, which would 
allow for a first-order phase transition for appropriate Higgs masses, we find \Z(ip,T)\ <; 1 for 
small if. We interpret this fact as an indication that higher loop-corrections as well as higher 
derivative terms in the effective action might be crucial for an understanding of the dynamics 
of the phase transition. For larger magnetic masses we obtain \Z(<p, T)\ ^ 1, and conclude that 
in this part of parameter space perturbation theory seems to be reliable. However, for such 
values of the magnetic gauge boson masses the phase transition is presumably of second order. 

Our negative one-loop result for Z(<p,T) for the range < 7 ^ 2 seems to indicate a 
decrease of the strength of the first-order phase transition. A further clarification of this point 
would however require the knowledge of higher loop-corrections to Z(ip,T) and perhaps also 
higher derivative terms in the effective action. 

Acknowledgment: We acknowledge valuable discussions with Scott Dodelson and Eric Wein- 
berg. 
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Appendix 



A The top quark self— mass 



In this appendix we evaluate the top quark self-mass S t . Among the various one-loop Feynman 
diagrams contributing to S t the dominant contribution comes from the gluon loops^. The 



corresponding contribution to Ej has been calculated in [|5| for the symmetric phase. For 
SU C (3) one obtains 

6 k 



t gluon 



(fco,0) 



-7o 



(a.i; 



Hence the gluon loops do not induce a plasma mass but an effective chemical potential n(ko,k = 
0) = g^-- The absence of a plasma mass in the symmetric phase is not surprising, because we 
do not expect chiral invariance to be broken to any finite order in perturbation theory. For the 
same reason the SUl(2) gauge boson loop-corrections to S 4 should not induce a plasma mass. 
Since furthermore g << g s , they will be neglected in the following, and we will concentrate on 
the scalar loop-corrections to 

The Higgs boson loop-contribution is given by 



Po 



d 3 p + m-t) 



(A.2) 



(27r) 3 p 2 -mf (p-ky ,,. r 

where po = (2n + l)iixT and ko = (2m + l)mT are fermionic Matsubara frequencies. The 
determination of a plasma mass or chemical potential requires the knowledge of £t(&o, 0): 



VjHiggs 



(fco,0) 



with 



h 



4 /t 



n=0 

00 

n=0 



70 / 2 

m< 2X" ^ ' 



<i 3 p 



(m t 2 - m 2 , + fc 2 J 



h 



8h 



■ftloh 



(A.3) 



(2vr) 3 (p(-) 2 -m?)(p(+) 2 -m 2 . 



rf 3 p / 1 

(2^)3 ^(-)2 _ m 2 



(A.4) 



(+)2 — fpp' 



Fermionic four-momenta are denoted by = {pi \p) whereas bosonic momenta are de- 
noted as p( + ) = (j)q,p) with Matsubara frequencies p$ = (2n + l)inT, p^ = 2nmT. The 
evaluation of I2 in the high-temperature expansion is straightforward, whereas I\ is somewhat 
more involved and given in appendix [FJ. Using these results we obtain to leading order in f t , 
Vb = m^/T and y f = m t /T 



^Higgs 



(fco,0) 



f-2 f 

Jt Cy f T+ h 



327T 



32tt 



To 



-Ck - k — 

2 7T Kq 



2 71 



(A.5) 



9 For a very heavy top quark with m t 
since then f t ~ g s . 



200 GeV one expects the Higgs loop to give a comparable contribution, 
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with C = -1 + 7b + In 2. 

The physical degrees of freedom are found by determining the poles of the propagator 
Sp{ko, k) = [| -mj- £i(&o, k)]~ l for k = 0. A straightforward calculation reveals four different 
zeros of Sp . In the symmetric phase, i.e., for yj = y b = 0, the pole positions are degenerate 
and given by ko = ±/fT/8, which is in perfect agreement with the findings of |^5, f|6|]0. 

The form of Et(&o,0) shows an additional plasma correction to the effective chemical po- 
tential as well as a plasma mass for the top quark of the form 



ft 



Vb 



ko 



n 



2 C "7 



kn 



5m t 



f, 



h Cy f T . 



(A.6) 



32vr V2 TV J k 32vr V2 tx J 32tt 
As expected, Sm t vanishes in the symmetric phase. Furthermore one easily sees that it is exactly 
canceled by the neutral Goldstone boson loop-contribution to E t , which may be calculated 
together with the charged Goldstone contributions analogously to Ef lggs . The final result is 



1 q 2 T 2 
6 k 



+ 



+ 



Jt 
32tt 

Jt 
32tt 



7o 



To 



1 T 2 / 1 

Ck y b k - — 1 + -y b 

71 k \ 71 

1 T 2 / 1 

2Ck y b k - — 2 + -y b 

TX K V TX 



1+75 

2 

1-75 



(A.7) 



B Some useful integrals 



In this appendix we provide the high-temperature expansion of the integral I± defined in ap- 
pendix |A] and for some related integrals also containing fermionic as well as bosonic momenta^]. 
Converting the Matsubara frequency summation into a contour integral and reading off the 
residues, one finds after evaluation of the spatial angular integral 



with Ay 5 



ut - Vf 

il b/f \y,z 



7T 



lj?\y b ,z b )-(Ay* + 7x 2 )l?>(y f 



2 and z 2 b/f = ±±Ay 2 + ^ (Ay 2 ) 2 + y 2 b and 



x 



n—l 



o + z 2 ^ y/x' 2 + y 2 (l =f exp ^fx r 



(b.i; 



(B.2) 



The integrals I^^(y, z) may be evaluated in the high-temperature expansion, i.e., for y, z << 
1, as follows: In a first step one determines l} 6 . Subsequently 1^1(1 * s calculated for y — 
and positive k. The leading y-dependence for arbitrary k may then be recovered by applying 



iS m (v, 



1 d 



r(b/f) 



o ± n+2 

ny ay 



(y, z) . 



(B.3) 



10 Actually the degeneracy already occurs for yf = but yb ^ 0. Furthermore it is lifted even in the symmetric 
phase for @. 



n Note that ' — k,jj ' is bosonic, though p 1 ^ ' and fcg ' are fermionic. 
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In this way one obtains 

7T 7T \ 47T/ 7T V 7T Z / y 

+ {-^+M- 2 -^- 4to £ + >)k 

T f&1 7T 1 7T 1 , ^ 2 

4 6; = 1 7s + n 2 + - 

3 2 2 ' 4 2 7r 



+ 



/F > = 7s - - In 2 y J 

4 2 2 



2 2 

4' = T ( 7 , + ]»2- T + -) 

,</> * 2 I 2 T 2 , „1 

+ i{i+7 E + 7 C -f-^ + ln2}, 2 + ?{ 7E -| + ln2}^, 
where 7^ ~ 0.577 denotes Euler's and 7 C ~ 0.916 Catalan's constant. 
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Figure 1: The qualitative shape of the effective potential is plotted for three different tempera- 
tures: T = T c , T b < T < T c and T = T b . Here T b denotes the barrier temperature for which the 
barrier between the two minima vanishes. For T b < T < T c we have also indicated the position 
of f mm and the potential difference AV(T) between the two minima. 



Figure 2: The total wave function correction factor Z(ip,T) (2. a), its scalar (2.b) and vector 
(2.c) field contribution are plotted for m_ff(T = 0) = 85 GeV, m t (T = 0) = 160 GeV and 
7 = 0, 1, 2 at the corresponding critical temperature T = T c . Note that T c is a function of '7, so 
that also Z sca i ar ((/?, T c ) is implicitly ^-dependent. 



Figure 3: The plotted boundary curve divides 7-^ space into two regions. Above the boundary 
\Z(ip,T)\ ~ \Z vcctOT ((p, T)\ < 1. Below it \Z(<p,T)\ > 1 and the wave function correction term 
dominates over the classical kinetic term of the effective action. 
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